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ABSTRACT 


Current  experiments  conducted  at  Dornier  System 
on  special  solid-surface  treatment  have  prompted,  a 
study  on  the  fluid  dynamic  implications  of  perfect- 
slip  flow.  Numerical  results  on  incompressible  fluid 
flows  past  thin  elliptic  cylinders  at  various  angles  of 
attack  and  at  moderate  Reynolds  numbers  show  that 
flow  separation,  instability,  and  vortex  shedding  occur 
also  under  the  perfect-slip  condition.  Of  particular 
interest  is  the  generation  and  spreading  of  an  initial 
vortex  after  the  abrupt  start  of  an  airfoil.  The  circu¬ 
lation  theory  of  lift  for  potential  flows  is  also  applicable 
to  perfect  slip.  Thus,  lift  exists  so  that  flying  under 
perfect  slip  is  possible. 


ADMINISTRATIVE  INFORMATION 

This  study  was  funded  by  the  Naval  Ship  Systems  Command  under 
the  Mathematical  Sciences  Program,  Subproject  SR  014  03  01,  Work 
Unit  1-1843-347. 

This  report  is  a  reprint  of  the  paper  "Laminar  Flow  Past  Bodies 
Undfer  Perfect  Slip,  "  published  in  the  proceedings  of  the  "Symposium 
liber  Gas-Oberflachen-Wechselwirkung,  "  May  3-5,  1972,  Meersburg, 
Lake  Constance,  Germany,  sponsored  by  the  German  Ministry  of  Defense 
and  organized  by  Dornier  System,  Friedrichshafen. 
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INTRODUCTION 


A  few  years  ago  Sydney  Goldstein  remarked  that,  if  a  fluid  could 

slip  freely  over  the  surface  of  a  solid  body,  the  world  would  be  very 

different.  *  Would  it  really  be?  To  give  a  partial  answer  to  this 

question  some  computer  experiments  have  been  made  which  simulate 

a  viscous  fluid  flowing  past  a  cylindrical  body  under  the  perfect-slip 

condition.  Perfect  slip  is  the  limiting  case  of  slippage  when  the  shear 

stress  at  the  solid  surface  vanishes.  This  is  a  hypothetical  situation 

since  slip  flow  of  a  Newtonian  fluid  still  awaits  verification  in  reality. 

Only  the  availability  of  electronic  computers  enables  us  to  simulate 

such  an  imaginary  flow.  It  is  not  the  purpose  of  this  paper  to  argue 

about  whether  and  how  slip  flow  can  be  obtained  within  the  realm  of 

continuum  physics  where  the  mean  free  path  is  very  small  with  respect 

to  the  body  length.  Apart  from  problems  in  the  rarefied-gas  regime 

or  from  certain  cases  of  non-Newtonian  liquids  (for  instance,  Bingham 

bodies)  where  slippage  is  known  to  exist,  a  renewed  interest  has  arisen 

in  how  to  obtain  slip  flow  for  a  Newtonian  continuum.  An  outgrowth  of 

this  effort  is  the  present  Symposium.  Interest  in  such  a  slip  flow 

was  also  the  incentive  for  this  study.  Some  numerical  results  on 

vorticity  production  and  propagation  under  perfect  slip  have  been 
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published  already. 


References  are  listed  on  page  17. 
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VORTICITY  PRODUCTION  ON  THE  BODY 

In  a  homogeneous  incompressible  fluid  without  nonconservative 
forces  acting,  vorticity  can  be  produced  only  at  the  body  surface.  By 
using  intrinsic  coordinates  s,  n  for  plane  or  axisymmetric  motions, 
where  s  coincides  with  the  streamlines,  and  n  is  orthogonal  to  s, 
the  vorticity  at  and  the  shear  stress  r  on  a  streamline  are  functions 
of  the  velocity  u  along  n  -  const: 

--  +  K  u  ,  (1) 

t/h  =  ku.  ,  (2) 

where  ju  is  the  dynamic  viscosity  and  K  the  curvature  of  the  streamline. 
Under  the  nonslip  (Condition,  u  s  0  at  the  body  surface  n  =  np 
Equations  (1)  and  (2)  reduce  to  co  =  r/ju  -  -  du/an.  Under  the  perfect- 
S lip  condition ,  r  ^  0  on  n  =  n^,  the  vorticity  is  co  =  2  K  u.  Hence, 
vorticity  occurs  at  the  surface  under  perfect  slip  only  for  k  ^  0.  Since 
every  cylindrical  body  with  a  closed  line  contour  of  its  cross-section, 
and  every  axisymmetric  body  with  a.  closed  surface  contour,  must  have 
areas  (or  points)  with  nonvanishing  curvature,  vorticity  will  always  be 
generated  under  perfect  slip  at  those  bodies. 

The  rate  of  transport  of  vorticity  from  the  surface  n  =  n^  into 
the  fluid  is  proportional  to  (aa>/dn)j.  This  term,  and  not  the  vorticity 
ujj  itself,  enters  the  equation  for  the  surface  pressure  p^  and,  hence, 
participates  in  determining  drag,  lift,  and  torque.  The  relation  between 
Pl  and  (dco/dn)i  follows  from  the  s -component  of  the  Navier-Stokes 
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equation: 


P1  pconst  "  £  £p^t'  +  u  Ss^  +  P  an  ^ds  *  ^ 

surface 

where  p  is  the  density  of  the  fluid  and  t’  the  time.  Under  perfect  slip 
the  term  (dto/dn)^  can  differ  from  zero  on  a  flat  part  of  the  surface, 
where  K  =  0  and  therefore  =  0.  In  this  case  vorticity  inside  the 
fluid  is  provided  from  other  areas  of  the  surface  where  K  /  0.  The 
only  case  for  which  coj  /  0  but  (dto/dn^  s  0  is  the  solid-body  rotation. 
No  exchange  of  vorticity  takes  place.  The  vorticity  is  ’’frozen”  inside 
the  fluid.  There  is  no  dissipation. 

Two  well-known  analytic  solutions  may  illustrate  the  role  of  vorticity 
production  in  determining  drag.  For  the  steady  flow  past  a  sphere  of 
unit  radius  at  vanishing  Reynolds  number  Re  =  2U p/p,  the  stream 
function  0  and  the  azimuthal  vorticity  w  are 

-  ~2  sin  0  (Ar  +  7  +  r  )>  w  =  — 2  sm  0  ,  (4) 

where  r,  p  are  the  spherical  polar  coordinates  in  the  meridional  plane. 

U  is  the  velocity  at  infinity.  The  integration  constants  A  and  B  are 
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A  =  -3/2,  B  =  1/2  for  nonslip  and  A  =  -1,  B  =  0  for  perfect  slip. 

The  drag  coefficient  (the  sum  of  C^F  and  C^p,  where  and 

CDp  are  the  portions  due  to  friction  and  pressure,  respectively)  is 

2 

defined  by  =  Drag/(p/2)7rU  .  Then,  for 

Nonslip:  C^Re  =  C^pRe  =  ^  ^DP^6  ~  ®  > 

Perfect  slip:  Cj-jRe  =  16,  C^Re  =  0,  C^pRe  =  16. 

Hence,  the  total  drag  coefficient  for  perfect  slip  is  tremendous,  2/3 
that  for  nonslip.  C^p  itself  is  even  twice  as  large  for  perfect  slip 
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although  (&co/dr)i  is  smaller.  The  larger  C^p- value  is  due  to  the 

contribution  of  the  surface  velocity  in  Equation  (3). 

If  the  sphere  is  flattened  to  an  oblate  spheroid,  an  analytic  solution 

3 

for  vanishing  Reynolds  number  exists  (Oberbeck  solution  ).  In  oblate 
spheroidal  coordinates  77 ,  0  for  the  meridional  plane  with 
x  +  iy  =  a  sinh  (77  +  iS),  ip  and  a)  are,  for  nonslip, 


U  £oslfe. 

2  cosh  77j 


K  v  .  2  - 
cosh  r? 


sinh  f^-1  ^ 

— — 2 -  cot  sinh  77)]. 

cosh  77^ 


2U  sin  0 
w  =  - - 

K  cosh  77  j 


tanh  n 

2  2  -  * 

cosh  rj-  sin  0 


(5) 


where  the  focal  distance  a  is  chosen  in  such  a  way  that  the  oblate  spheroid 
has  unit  radius,  that  is,  a  =  1/cosh .  The  surface  of  the  body  is 
r)  =  fjj  .  The  constant  K  is 


1  -  2  -  -1  - 
K  =  - 2 —  [  sinh  Vi  ~  (sinh  77^  -  1)  cot  sinh  77^  ] . 

cosn  77j 

For  the  flat  disk  7^  =  0  we  observe  the  following  curious  situation: 
Except  at  the  edge  the  surface  vorticity  is  zero  and,  therefore,  the 
shear  stress  is  zero,  whereas  (dco/dT?)^  /  0  in  general.  The  total 
vorticity  field  is  generated  at  the  focal  point  which  for  77^  =  0  coincides 
with  the  edge.  The  situation  is  illustrated  in  Figure  1  for  the  fat  oblate 
spheroid  77^  =  0.5  to  show  the  lines  of  constant  vorticity  inside  the  body. 
Notice  that  (Sco/Sr? ) ^  changes  its  sign  along  the  surface.  Near  the 
centerline  vorticity  is  transported  from  the  fluid  to  the  body.  Around 
the  edge  vorticity  emanates  from  the  surface.  For  77^  =  0  this  region 
is  confined  to  the  edge. 
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Figure  1.  Lines  of  Constant  Vorticity  for  an  Oblate  Spheroid 
rj  =  0.  5  at  Re  =  0  Under  Nonslip 

Under  the  perfect-slip  condition  the  solution  is  not  separable,  and 
an  analytic  integral  could  not  be  obtained.  A  numerical  solution  was 
therefore  constructed  for  =  0.  05.  (The  numerical  procedure  breaks 
down  for  =  0. )  The  drag  coefficients  are: 

Nonslip,  =  0  :  C^Re  =  20.  372,  C^p  =  0, 

Nonslip,  77^  =  0.  05:  C^Re  =  20.  395,  C^p  =  1.  504, 

Perfect  slip,  77.^  =  0.05:  CpRe«28  ,  C^p  *  0. 

Although  the  numerical  value  for  perfect  slip  is  very  crude,  the  result 

shows  that,  as  in  the  case  of  the  sphere  flow,  C^p  is  greater  for 
perfect  slip  than  for  nonslip. 


FLOW  BEHAVIOR  UNDER  PERFECT  SLIP 

The  study  of  flows  past  bodies  for  nonvanishing  Reynolds  numbers 

requires  the  use  both  of  numerical  methods  and  of  electronic  computers. 
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An  existing  computer  code  for  moderate  Reynolds. numbers  has  been 

employed  to  calculate  two-dimensional  laminar  flows  of  an  incompressible 

fluid  past  elliptic  cylinders  under  perfect  slip.  Some  data  have  been 
2 

published  already  so  that  the  results  reported  here  may  be  considered 
as  complements. 

In  this  computer  code  the  elliptic  cylinder  is  imbedded  in  a  network 
of  points  which  are  the  intersections  of  constant  elliptic  coordinate  lines. 
The  elliptic  coordinates  77, 0  are  related  to  the  Cartesian  coordinates 
x,y  by 

x  +  iy  =  a  cosh (77  +  i0  ),  a  >  0,  (6) 

where  a  is  the  focal  distance.  The  body  contour  is  given  by  77  =  77^ 
(Figure  2). 


Figure  2.  Elliptic  Coordinate  System  and  Definition  of  Angle 

of  Attack  a  . 
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The  Reynolds  number  is  defined  by  Re  =  dUp/ji  with  d  =  a  cosh  . 

The  calculation  starts  with  a  potential  flow  as  initial  condition  at 
t  =  t’U/a  =  0„ 

The  following  examples  have  been  selected: 

0!  =  0°,  Re  =  200,  17^  =  0. 1,  symmetric,  almost  steady  flow 
at  =  90°,  Re  =  10,  172  =  0- 1,  symmetric,  almost  steady  flow 
a  ~  45°,  Re  =  200,  77^  =  0. 05,  asymmetric, transient  flow 

Flows  past  thin  elliptic  cylinders  under  perfect  slip  show  the  largest 
divergence  from  those  under  nonslip  when  ot  =  0°.  The  reason  is  that, 
under  nonslip,  friction  drag  predominates  over  pressure  drag.  Under 
perfect  slip,  however,  friction  drag  is  zero  by  definition.  For  a  =  0°, 
Re  =  200,  77^  =  0. 1  the  drag  coefficient  under  perfect  slip  is  1/60  of  that 
under  nonslip  ,  and  the  surface  velocity  v.  (made  dimensionless  by  U) 
comes  quite  close  to  that  of  potential  flow  (see  Figure  3).  The  velocity 
profile  normal  to  the  surface  at  0  =  7r/4,  where  the  77 -coordinate 
coincides  with  the  y-axis  in  Figure  2,  reveals  that  flow  under  perfect  slip 
approaches  the  potential  flow  away  from  the  body  much  sooner  than  flow 
under  nonslip  (Figure  4).  This  indicates,  if  extrapolation  to  larger 
Reynolds  number  is  allowed,  that  the  boundary  layer  under  perfect  slip 
is  thinner  than  that  under  nonslip. 

Blunt  bodies  or  plates  perpendicular  to  the  main  flow  behave 
differently.  As  already  indicated  in  the  solutions  for  Re  =  0,  the  drag 
coefficients  for  motions  under  perfect  slip  and  under  nonslip  are  of  the 
same  order.  Figure  5  shows  the  transient  stage  of  a  flow  normal  to  a 
thin  elliptic  cylinder.  The  drag  coefficient,  which  is  defined  here  by 
=  Drag/(p/2)dU^,  is  plotted  against  t  for  a=  90°,  Re  =  10,  772  =  0. 1. 
Only  a  slight  difference  between  nonslip  and  perfect  slip  is  observed. 

The  surface  velocity,  which  is  plotted  in  Figure  6,  deviates  considerably 
from  that  of  potential  flow.  The  occurrence  of  a  negative  velocity  in  the 
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Figure  3.  Surface  Velocity  for  a=  0°,  Re  =  200,  n ^ 
at  Almost  Steady  State 


Figure  4.  Velocity  Profiles  for  a  =  0  ,  Re  =  200 
0=  7t/4  at  Almost  Steady  State 


Figure  5.  Drag  Coefficient  Against  Time  for  a  =  90  ,  Re  =  10, 

Abrupt  Start  at  t  =  0. 


wake  points  to  flow  separation,  which  is  verified  in  Figure  7.  Actually, 

the  flow  patterns  in  Figure  7  are  almost  the  same  as  those  for  motions 
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under  nonslip  (see  streamline  and  vorticity  patterns  )i  This 

result  is  surprising  since  we  expect  perfect-slip  flow  to  behave  much 

like  potential  flow.  The  development  of  the  wake  is  recorded  in  Figure  8, 

where  the  dimensionless  wake  length  s’/d  is  plotted  against  time. 

6  7 

Experimental  data  from  Taneda  and  Taneda  and  Honji  are  included  in 
the  graph.  Their  results  are  discussed  by  Liigt  and  Haussling.  ’ 


GENERATION  OF  LIFT  UNDER  PERFECT  SLIP 

Classical  aerodynamics  explains  the  generation  and  maintenance  of 
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lift  of  an  airfoil  by  means  of  the  Lanchester-Prandtl  hypothesis  :  After 

the  abrupt  start  of  the  body  a  vortex  at  the  trailing  edge  develops  and 

separates,  leaving  behind  a  bound  vortex  of  equal  strength  but  opposite 

rotation  around  the  airfoil.  Mathematically,  a  circulation  around  the 

airfoil  is  superposed  on  the  potential-flow  solution.  This  circulation 

gives  rise  to  a  lift  force  whose  magnitude  is  determined  by  the  Kutta 

condition  at  the  trailing  edge.  Since  no  consideration  of  viscosity  is 

required  in  this  model  (in  contrast  to  the  drag  force),  the  nonslip  condition 

of  a  viscous  fluid  cannot  be  applied.  The  model  has  been  proved  to  be 

remarkably  good  for  high  Reynolds -number  flows.  Even  for  Reynolds 

numbers  as  low  as  200  the  essential  features  of  the  model  are  observed.  ^ 

If  it  can  be  shown  that  flow  under  perfect  slip  satisfies  the  Kutta 

condition  (in  the  form  extended  to  viscous  fluids),  the  Lanchester-Prandtl 

hypothesis  must  also  be  valid  for  perfect-slip  flow.  The  fulfillment  of 

2 

the  Kutta  condition  for  perfect-slip  flow  has  been  demonstrated 
for  Re  =  200,  a  -  45°.  Assuming  that  this  result  is  valid  for  higher 
Reynolds  numbers  too,  we  arrive  at  the  important  conclusion  that  the 
classical  airfoil  theory  is  applicable  to  perfect  slip  flow.  This  statement 
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covers  three-dimensional  effects  since  the  trailing  vortices  at  the  tips 
of  the  wings  exist  also  under  perfect  slip.  Thus,  an  airplane  can  fly 
under  perfect  slip! 

Let  us  consider  the  initial  phase  of  the  process  during  which  the 
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Kutta  condition  is  established  and  which  is  not  covered  by  Lugt.  The 

following  description  holds  for  both  nonslip  and  perfect  slip. 

At  t  =  0  a  potential  flow  around  a  thin  elliptic  cylinder  for  a  =  45°  is 

assumed  without  circulation  (Figure  9),  There  is  no  lift  (and  no  drag). 

At  this  instant  the  initial  irrotational  disturbance  caused  by  the  sudden 

motion  of  the  body  is  felt  immediately  everywhere  but  vorticity  has  not  yet 

started  to  propagate.  The  generation  of  a  discontinuity  sheet  for  the 

vorticity  on  the  body  surface  at  t  =  0  causes  an  infinite  lift  and  drag  force. 

Subsequently,  the  sheet  spreads  out  to  a  boundary  layer,  and  the  lift  and 
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drag  force  drop  to  finite  values.  (In  Lugt,  -  'last  paragraph,  a  difference 
was  reported  between  perfect  slip  and  nonslip.  It  was  found  that  the 
numerical  result  for  perfect  slip  just  after  t  -  0  is  not  physically  realistic 
but  is  the  result  of  the  numerical  scheme.) 

With  time  inertia  carries  vorticity  fromithe  boundary  layer  down 
stream.  This  causes  the  rear  separation  point  of  the  zero  streamline 
to  move  towards  the  trailing  edge.  Thus,  the  Kutta  condition  becomes 
fulfilled  at  about  t  =  0.  4  (Figure  9).  For  sufficiently  high  Reynolds 
number  (observed  for  Re  =  200)  the  vorticity  field  forms  a  local  extremum 
behind  the  trailing  edge  which  may  be  interpreted  as  the  starting  vortex 
(Figure  10).  However,  for  Re  =  15  and  30  under  nonslip,  no  extremum 

4 

of  the  vorticity  is  observed,  although  the  Kutta  condition  is  satisfied. 
Diffusion  of  vorticity  dominates  here  over  convection.  The  Lanchester- 
Prandtl  hypothesis  does  not  hold  for  this  low  Reynolds -number  flow.  It 
must  be  emphasized  that  for  all  Reynolds  numbers  the  ultimate  cause  of 

lift  (and  drag)  is  the  generation  of  the  vorticity  sheet  at  t  =  0. 

The  streamlines  associated  with  the  vorticity  extremum  are  wavy 
(Figure  9).  They  represent  a  vortex  moving  relative  to  the  reference 
frame,  that  is,  relative  to  the  body. 
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PERFECT  SLIP 


t=0 

Hi=0 


Figure  9.  Sequence  of  Streamline  Patterns  for  a  *  45°,  Re  =  200  at  Various 
Times  Immediately  after  the  Abrupt  Start  at  t  =  0. 


CONCLUSIONS 


Perfect-slip  flow,  if  it  can  ever  be  approached  in  reality,  differs 
drastically  from  nonslip  flow  when  the  motion  under  the  nonslip  condition 
is  such  that  the  shear  stress  on  the  solid  wall  dominates  over  the  pressure 
force.  Slender  bodies  parallel  to  the  main  stream,  fluid  motion  in 
capillaries  and  pipes,  and  rotating  devices  such  as  ball  bearings,  would 
benefit  most  from  the  slip  condition  in  that  drag  and  pressure  loss  would 
be  greatly  reduced.  On  the  other  hand,  flows  past  blunt  bodies  (and 
large  displacements  of  fluid  in  general)  are  not  greatly  affected  by  the 
change  in  the  surface  condition.  Classical  airfoil  theory  cannot  distinguish 
between  nonslip  and  perfect  slip  since  the  Kutta  condition  is  satisfied  for 
both  types  of  motions.  Flow  separation,  vortex  shedding,  and  instability 
(which  is  not  explicitly  mentioned  in  this  paper,  but  which  precedes 
vortex  shedding)  are  phenomena  which  also  occur  under  perfect  slip. 

This  implies  that  turbulent  boundary  layers  at  high  Reynolds  numbers 
must  exist.  Turbulence  generated  away  from  solid  surfaces  is,  of  course, 
not  affected  by  the  slip  condition.  Turbulence  dominates  in  fluid  phenomena 
in  technology  and  in  our  every  day  life.  Thus,  returning  to  Goldstein’s 
remark  cited  in  the  introduction,  that  the  world  would  be  very  different 
under  slippage,  perhaps  it  would  be  to  a  certain  extent,  but  not  as  much 
as  Goldstein  might  have  envisioned.  Until  perfect-slip  flow  is  verified 
experimentally,  its  implications  for  the  world  in  which  we  live  would 
make  a  good  subject  for  science  fiction. 
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